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Abstract. We present our preliminary results on the asteroseismological investigations of the ZZ 
Ceti star GD 244. We used literature values of the effective temperature and surface gravity and 
utilized the White Dwarf Evolution Code of Bischoff-Kim, Montgomery and Winget (2008, ApJ, 
675, 1512) to build our model grid for the seismological analysis. Five observed pulsational modes 
published up to now were used to find acceptable model solutions. We found that the best model 
fits have masses between 0.61 and O.74M 0 and constitute two groups with hydrogen layer masses 
of either ~ 10 5 or 10 6 M . : . Based on a statistical analysis of a larger sample of possible model 
solutions, we assume that the mass of the star is below ~ 0.68 M 0 and the oxygen content in the 
centre is less than 60 per cent. 
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INITIAL PARAMETERS AND THE MODEL GRID 

Our aim is to find stellar models that match the observed properties of the ZZ Ceti star 
GD 244 with high precision. We present our asteroseismological results using published 
period and atmospheric parameters. 

GD 244 was observed with the Canada-France-Hawaii Telescope (CFHT) in 1999 
(1 night, [3]) and at McDonald Observatory in 2003 (10 nights, [8]). Table 1 shows the 
period and amplitude values detected and the assumption about mode identification. Not 
all the periods were observed in both data sets. Amplitude variations may be responsible 
for the differences. In such cases we could attempt to find normal modes enough for 
asteroseismology by taking all known modes from the different seasons (see e.g. [5]). 
Period values used for the seismological analysis are typeset in boldface. All but one 
were determined from the McDonald Observatory data set, which had better frequency 
resolution than the short CFHT data. In three cases the same m = -1 value was presumed 
for the modes. 

Table 2 summarizes the atmospheric parameters determined by optical spectroscopy. 
We covered the log g range and a wider range in T e ff with our models. The latter 
approximately corresponds to the interval of ZZ Ceti stars. 

We used the White Dwarf Evolution code of Bischoff-Kim, Montgomery and Winget 
[1] to build our model grid. We varied five input parameters of the WDEC in the 
following ranges (in square brackets: step sizes): 


TABLE 1. Pulsation modes of GD 244. A [No.] refers to 
relative values estimated by Fig. 3 in [3], 
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TABLE 2. Effective temperature and surface 
gravity values determined by spectroscopy. 


Teff [K] 

logg [cgs] 
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T eff = 10 800 - 12 200 [200] K, 

M* = 0.525 - 0.740 [0.005] M 0 (log g ~ 7.85 - 8.23), 

M n = 10~ 4 - 1(T 8 [10' a2 ] M,, 

Xo = 0.5-0.9 [0.1] (central oxygen abundance), 

Xf m = 0.1 - 0.5 [0.1] (the fractional mass point where the oxygen abundance starts 
dropping). 

We fixed the mass of the helium layer at 10 2 A/.,. 

RESULTS AND DISCUSSION 

Our first criterion for an acceptable model was that it should give period values close to 
the observed ones. We used the parameter root main square ( r.m.s .) calculated from the 
observed and model periods to select among the possible solutions. The fitting routine 
FITPER [4] was applied for this purpose. Since we had only five modes to fit, many 
acceptable models had low r.m.s. value. We let all five modes to be l — 1 or 2 for the 
fitting procedure. Assuming better visibility of / = 1 modes, we selected the models that 
gave at least three l — 1 solutions. Beside the low r.m.s., this was our second criterion. 

Table 3 summarizes the parameters of our best-fitting models. The last column shows 
the 1 — 2 modes only. We can discriminate between the models based on the observed 
amplitudes. Considering that the 256 and 307-s modes have the largest amplitudes, it 
is improbable that both of them are l — 2. Thus, the 0.665 and 0.685 A/ ., models (both 
at the low temperature region) are less favourable. The 0.620 M Q model represents an 
interesting case. It has the thinnest hydrogen layer and only this gives l — 1 value for 
both the dominant modes. The rest of the models form two groups: one with 10~ 6 M* 



TABLE 3. Parameters of our best-fitting models. We indicate the 
r.m.s. values calculated from the observed and model periods in 
column 6. The last column shows the modes of Table 1 that are / = 
2 according to the model. The other modes are 1 = 1. 
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hydrogen layers and 1—2 values for the 256 and 294-s modes, and another with 
~ 10 5 M* hydrogen layers and 1 — 2 value for the 307-s mode. Since the 307-s mode has 
the largest light amplitude, being / = 2 would imply a much larger physical amplitude, 
thus we prefer an / — 1 solution for this mode. Therefore, the models with stellar masses 
between 0.610 - 0.630 M & and Mu — 10or the one with Mu = 1.6 x 10~ 7 M* are 
better choices. The 203-s mode is always l — 1. 

Yeates et al. [8] also identified the 203-s mode as l — 1 using combination frequency 
amplitudes. Castanheira & Kepler [2] gave 1 — 2 value for this mode. They built their 
model grid using fixed, homogeneous C/O 50:50 cores, but varied the mass of the helium 
layer. Only the 203, 256 and 307-s modes were fitted. Their best solution was T e ff = 
12 200 K, M* - 0.68 M Q ,Mu = 10” 7 M* and M Us = 10~ 3 ' 5 M*. 

We investigated the acceptable models from a statistical point of view also. We 
extended our sample selecting every model that gave at least three l — 1 solutions and 
had r.m.s. <1.5 values. With these criteria we obtained significantly more, 81 models 
for the analysis. The histograms in Fig. 1 show the number of models in the bins given 
by the step sizes for the five physical parameters. 

As Fig. la shows, most of the solutions are in the 11 400 - 12 000 K T e ff range and 
the most frequent value is 11 600 K. Most of the spectroscopic T e ff values are also 
around 11 600 K. Considering Fig. lb, the populated bins can be found between 0.62 and 
0.67 M & . Peak values are at 0.63 and 0.65 M Q . Fig. lc clearly shows the two families of 
solutions with Mu ~ 10 5 and 10 6 AT*. The favoured hydrogen layer mass with this 
method is 10~ 5 M*. According to Fig. Id, the central oxygen abundance could be 50 - 
60% (or less), larger values are not preferred. We can not give such constraint on the Xf m 
parameter (Fig. le). 

Additional photometric observations were obtained on this star both at McDonald 
Observatory and Piszkesteto mountain station of Konkoly Observatory. The analysis 
of these data and a detailed asteroseismological investigation will be the subject of a 
forthcoming paper. 
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FIGURE 1. Histograms of the 81 selected models for the five physical parameters varied in the grid. 
The bin sizes are suited to the grid step sizes. 
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